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Abstract Thermohaline mixing has been recently identified as the probable domi- 
natin g process that governs the p hotospheric composition of low-mass bright giant 



stars jCharbonnel & Zahnl 120071) . Here we present the predictions of stellar mod- 
els computed with the code STAREVOL including this mechanism together with 
rotational mixing. We compare our theoretical predictions with recent observations. 



1 Introduction 

The standard theory of stellar evolution predicts that the convective envelope of low- 
mass stars deepens in mass during the contraction of the degenerate He-core after 
the main sequence turnoff, and engulfes hydrogen-processed material (the so-called 
first dredge-up, hereafter Idup). This induces a decrease of the surface '^C/'^C ratio 
and of the Li and '"C abundances, while '^N and ^He abundances increase. After 
the Idup, the convective envelope retracts and the hydrogen burning shell (HBS) 
moves outward (in mass). According to the standard theory, no further change of 
the surface chemical properties is expected a fter the Idup on the Red Giant Branch 
(RGB). However, spectroscopic observations dGilrov , ll989HGilrov & Browni[l99T 



C. Charbonnel 

Observatoire de Geneve, Universite de Geneve, 51 ch. des Maillettes 1290 Versoix (Switzerland), 
Laboratoire d'Astrophysique de Toulouse-Tarbes, CNRS UMR 5572, Universite de Toulouse, 14 
av. E.Belin F-31400 Toulouse, 
e-mail : jchabonnel . CorinneBunige . ch| 

N. Lagarde 

Observatoire de Geneve, Universite de Geneve, 51 ch. des Maillettes 1290 Versoix (Switzerland), 
e-mail: [Lagarde .NadegeBunige . ch| 

P. Eggenberger. 

Observatoire de Geneve, Universite de Geneve, 51 ch. des Maillettes 1290 Versoix (Switzerland), 
e-mail : |Patrick . Eggenberger@unige . ch| 



1 



2 Corinne Charbonnel, Nadege Lagarde & Patrick Eggenberger 



Luckl 1994 


Gratton et al.l 2000|; Tautvaisiene et al., 2000|; Tautvaisiene et al.. 


2005; 


Smith et al. 


, 20021 Shetrond, 2003; Pilachowski et al., 2003; Geisler et all 


2005; 


Spite et al., 20061 Recio-Blanco & de Lavernvl 20071 ISmilianic et al., 2009) show 



clear signatures of "extra-mixing" on the upper RGB in low-mass stars, when the 
HBS crosses the discontinuity left behind by the Idup at the bump luminosity. More 
specifically, the carbon isotopic ratio and the C and Li abundances drop again as 
the stars move across the bump. Different processes have been proposed to explain 
these abundance anomalies. Here we recall the potential impact of rotation-induced 
mixing and thermohaline instability on the RGB. 



2 Rotation-induced mixing 



Rotation-induced mixing has an impact on the internal abundance profiles of chemi- 

cals involved in hydrogen-burning while the stars are on the main sequence dTalon & Charbonneli 
1998ll2010l;ICha rbonn el & Talonlll999l;|Palacios et al l 120031; IPasquini et al.l booi 
Smiljanic et alll2009i; iCharbonnel & Lagardei l2010l) . This can be seen in Fig. [T] 
where we present the chemical structure at the end of the main sequence for 1 .25 Mq 
stellar models computed assuming different initial rotation velocities. In the rotat- 
ing models, the abundance gradients are smoothed out compared to the standard 
(i.e., non-rotating) case; ^He, '^C, '^N, and '^O diffuse outwards, while '"C and 
'^O diffuse inwards. After the turnoff , this reflects in different post-dredge up 
predictions for the surface abundances that agree well with observations in low- 
luminosity red giant stars (see Fig. 2). However, rotation-induced mixing does not 
explai n the abundance anomahes observed in low- mass red giants brighter than the 
bump dChaname et al.L 120051; iPalacios et al. [ |2006h . 



3 Thermohaline mixing 



Charbonnel & ZahnI ( l2007h identified thermohaline mixing as the process that gov- 
erns the photospheric composition of low-mass stars above the bump on the RGB. 
At this evolution ary phase, this double d iffusive instability is induced by the ^He(^He; 
2p)^He reaction ( Eggleton et al.l l2006b that creates an inversion of mean molecular 
weig ht jUkichl l 197 lh. Here we use the following prescrip tion for the turbulent dif- 



fusivity coefficient dUlrichi 119721; iKippenhahn et al. ■ Il980l) : 



D, 



'ad - 



-V) 



forV^ <0, 



(1) 



Ct = -TT^a^, 



(2) 
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Fig. 1 Chemical structure at the 
tumoff of a 1.25 M j star com- 
puted for different initial rota- 
tion velocities as indicated. The 
mass fractions are multiplied by 
100 for ^He, '^c, and "*N, by 
2500 for "C, by 50, 900, and 
5x10"* for "'O, "O, and ''^O 
respectively, and by 1500 for 
^^Na. The vertical arrows show, 
in all cases, the maximum depth 
reached by the convective en- 
velope at its maximum extent 
during the Idup. Figure fro m 
ICharbonnel & Lagardel j2010l) . 



with K the thermal diffusivity and a = 5 the aspect ratio of sah fingers dUkich , 



1972h . At the RGB bump and above, the thermohaline diffusion coefficient is several 



order of magnitudes higher than the total diffusion coefficient related to rotation- 
induced processes (iCharbonnel & Lagai"deti201Q) . 



4 Model predictions and comparisons to observations. 



The re sults shown in this section are presented and discussed in detail in lCharbonnel & Lagarde 



(l2010l) and have been computed with the code STAREVOL including rotation 



induced processes and thermohaline mixing. Here we briefly discuss the effects of 
rotation-induced mixing and thermohaline mixing on carbon isotopic ratio, lithium 
and ^He. 



4.1 Carbon isotopic ratio 



In Fig|2]we compare the evolution of the theoretic al surface carbon isotopi c ratio for 
1 .25 Mq-j models with observations in M67 stars bv lGilrov & Brownl(ll99ll) . We note 
that rotation-induced mixing on the main sequence slightly lowers the post-dredge- 
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up '^C/'^C value (see §2). On the other hand, thermohaline mixing leads to further 
decrease of the carbon isotopic ratio at the luminosity of the bump (log (L/L, 7,) ^2), 
in excellent agreement with M67 data. 



Fig. 2 Evolution of the surface 
'^C/'^^C ratio as a function of 
stellar luminosity for models 
of 1.25 Mq stars. Different 
tracks are for different initial 
rotation velocities (non-rotating 
case, 50, 80, and 110 km.s"' 
respectively in black, red, 
green , and blue). Obs er vations 
by iGilrov & Browi] l ll99lh 
in evolved stars of the open 
cluster M67 (tumoff mass 
~1.2 Mi.) are also shown 
(triangle, squares, and circles 
for subgiant, RGB, and clump 
stars respectively). Figure from 
.Charbonnel & Lagarde C2010.) . 

In Figl3]we show the predictions for the surface ^^C/^^C ratio at the tip of the 
RGB and at the end of second dredge-up (black and blue lines respectively) for 
models over a large mass range and compare them with observations in stars be- 
longing to open clusters of various turnoff masses. We see that in stars with masses 
below '-^2M,3, thermohaline mixing is the main physical process governing the pho- 
tospheric composition of evolved giants, although the final carbon isotopic ratio also 
slightly depends on rotation-induced mixing on the main sequence. In intermediate- 
mass stars that do not reach the bump on the RGB and do not undergo thermoha- 
line mixing at that phase, rotation is necessary to explain the data and accounts for 
star-to-star abundance variations at a given evolutionary status. Overall, the present 
models explain very well the observed abundance patterns over the considered mass 
range. 




4.2 Lithium 

In FiglHwe present lithium data for field red giant stars with metallicities around 
solar and precise determination of their evolutionary status, and compare them to 
predictions for models of various masses. Contrary to the standa rd case, rotation- 



induced mixing leads to Li depletion on the main sequence (see e.g. lTalon & Charbonnel , 
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Fig. 3 Observ a tions o f '^C/'-^C in evolved s tars o f Galactic open cluste rs by ISmilianic et alj 
j2009h . iGilrovl jl989[) , lOilrov & Browr] il99lh . and iMikolaitis et alJ <2O10l) as a function of the 
tumoff mass of the corresponding host cluster. Squares, circles, and as teriscs are for RGB, clump, 
and early-AGB stars respectively, while diamonds are for stars from iGilrovl <1989l) sample with 
doubtful evolutionary status; triangles are for lower limits. A typical error bar is indicated. Theo- 
retical predictions are shown at the tip of the RGB and after completion of the second dredge-up 
(black and blue lines respectively). Standard models (no thennohaline nor rotation-induced mix- 
ing) are shown as dotted lines, models with thennohaline mixing only (Vzams=0) as solid lines, 
and models with thermohaline and rotation-induced mixing for differen t initial rotation veloci- 
ties as indicated as long-dashed, dot-dashed, and dashed lines. Figure from lCharbonnel & Lagard3 



19981120101: iPalacios et al.L 120031) . After the Idup, the surface Li abundance remains 



constant until the stars reach the bump luminosity where thermohahne mixing be- 
comes efficient and destroys Li. Later on the star reaches the RGB tip, and finally the 
second dredge-up decreases again Li at the surface. Models are in perfect agreement 
with observations. 

We have computed a few models along the TP-AGB, and found that thermohaline 
mixing leads to non negligible fresh lithium production, as shown in Fig|5] There 
we present the evolution of the surface lithium abundance as a function of both 
effective temperature and bolometric magnitude for TP-AGB models of 1 .25 and 
2.OM0 stars. Theoretical predictions are compared with observations of the sample 
of low-mass oxygen-rich AGB variables belonging to the Galactic disk studied by 



Uttenthaler & Lebzelter and are found to fit very nicely the observed Li 
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Fig. 4 Lithium data for field 
evolved stars from the sample by 
Charbonnel et al. (in prep.) that 
are segregated according to their 
mass (left and right panels in- 
clude respectively sample stars 
with masses lower and higher 
than 2 Mj,; Li detections and 
upper limits are shown as cir- 
cles and triangles respectively). 
Theoretical lithium evolution is 
shown from the ZAMS up to 
the end of the early- AGB. Var- 
ious lines correspond to predic- 
tions for stellar models of dif- 
ferent masses computed without 
or with rotation as indicated, and 
with thermohaline mixing in all 
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behaviour. Let us note that despite the strong production of fresh Li at that phase, 
the total stellar yields remain negative for this element. 




1.35Mg V2^^g=0km.s-' 
l.SSMj, V,„5=110km.s-i 

V„„=110km.s-' 




4000 3500 3000 2500 
Teff 



-3 -4 
Mbol 



Fig. 5 Lithium observations in 
oxygen-rich variables belong- 
ing to the Galactic disk as 
a function of effective tem- 
peratur e and bolometric mag- 
nitud e iUttenthaler & Lebzelteri 
I2OIOI) : circles and triangles are 
for abundance determinations 
and upper limits respectively. 
Theoretical lithium evolution is 
shown from the early-AGB up 
to the end of the TP- AGB. Var- 
ious lines correspond to predic- 
tions for stellar models of dif- 
ferent masses computed with- 
out or with rotation as indi- 
cated, and with thermohaline 
mixing in all cases. Figure fro m 
ICharbonnel & LagardeH2010l) . 
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4.3 Helium 3 

On th e main sequ ence, a ^He peak builds up due to pp-reactions inside low-mass 
stars ( IbenL 1967 ). and is engulfed in the stellar envelope during the Idup. As a 
consequence the surface abundance of ^He strongly increases on the lower RGB 
as can be seen in Fig|6] which presents the evolution of ^He mass fraction at the 
surface of 1 Mf, model at solar metallicity in the standard case and in the case 
with thermohaline mixing (black solid and red dotted lines respectively). After the 
bump, thermohaline mixing transports ^He from the convective envelope down to 
the hydrogen-burning shell where it burns. This leads to a rapid decrease of the 
surface abundance (and thus of the corresponding yield) of this element as can be 
seen in Fig|6] We are presently computing similar stellar models over a large range 
in stellar mass and metallicity in order to quantify the actual contribution of low- 
mass stars to Galactic ^He (Lagarde et al., in preparation). We are confident that the 
corresponding ^He yields will help reconciling t he primordial nucle osynthesis with 
measurements of ^He/H in Galactic HII regions ( CharbonnelL 2002 ). 



1.0z014_stondo : May 10 11:19 2010 
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Fig. 6 Evolution of tlie sur- 
face abundance of ^He (in 
mass fraction) from the 
pre-main sequence up to 
tlie AGB tip for IM.t, mod- 
els at solar metallicity. The 
black solid line and the 
red dotted-line correspond 
to the standard and thermo- 
haline cases respectively. 



Page: 7 job : Charbonneletal maoro : svmult . els date/time : 2-Sep-2011/0 : 14 



8 

5 Conclusions 
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An inversion of molecular weight created by the ^He(^He; 2p)^He reaction is at the 
origin of thermohaline instability in low-mass RGB stars brighter than the bump. 
The associated mixing explains very well observations of '^C/'^C and Li in these 
objects. Rotation-induced mixing, coupled with thermohaline mixing, allows us to 
explain the '^C/'^C and Li data in giant stars over a large mass range. Thermohaline 
mixing has also an effect during the TP-AGB phase, where it leads to fresh lithium 
production. Finally, this process helps reconciling stellar yields predictions with the 
Galactic evolution of ^He as constrained by the data in Galactic HII regions. 
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